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1 . Summary 


This semiannual report briefly summarizes research activities 
performed during the period of June 26, 1993 through August 31, 
1993. We developed the Robust Stability of Systems where transfer 
function or Characteristic polynomial are multilinear affine 
functions of parameters of interest in two directions, Algorithmic 
and Theoretical. In the algorithmic direction, a new approach that 
reduces the computational burden of checking the robust stability 
of the system with multilinear uncertainty is found. This 
technique is called "Stability by linear process." In fact, the 
"Stability by linear process" we describe here gives an algorithm. 
But we still have something else to be done. In analysis, we 
obtained a robustness criterion for the family of polynomials with 
coefficients of multilinear affine function in the coefficient 
space and obtained the result for the robust stability of diamond 
families of polynomials with complex coefficients also. We 
obtained the limited results for SPR design and we provide a 
framework for solving ACS. Finally, copies of the outline of our 
results are provided in the appendix. Also, there is 
administration issue in the appendix. 


2.1 Parametric Robust Stability. 

Let p (s, Q) = E" i=0 a i (q)S i , q e Q C R n 
and aj(q) are multilinear affine functions. As we know 
a[p(s,Q)] C C\ ** p ( s , Q) e H 0 £ p(ja>, Q) 



In this section, we develop 

p (s , Q) e H 4* 0 £ p(jw, Q) 
by algorithm, and 

o [p (s, Q) ] C C » p(s,Q) e H 
by theoretical analysis. 
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2.1.1 Stability by linear process. 

Our problem is to find a quick way to determine if the origin 
0 in the complex plane C within the image of hype rrec tangle Q C R n 
under the mapping 

p(jco,q) = f (q) + g(q)j for fixed co 
Where f (q) and g(q) are multilinear affine functions for V q e Q. 

We completely utilize linearity of multilinear affine function 
to study the line segments parallel to axis in Q mapping to the 
line segments in C. In accordance with our investigation, these 
line segments in C can be classified two kinds, one is so-called 
"perpendicular" and the other is so-called "non-perpendicular." 

Now we define that a line segment 1' in C is said to be a 
perpendicular if there is a point A' on 1' such that the line 
through 0 and A' , OA' perpendicular to 1' . If not, then this line 
segment 1' is called non-perpendicular. (see figure 1) 



Figure 1 

In Figure 1 (a) and (b) are the case of perpendicular; (c) is 
the case of non-perpendicular. 

It is evident that (1) if all 1' in C are perpendicular then 
0 e p ( j co , Q) , (2) if all 1' in C are non-perpendicular then 0 £ 

p(jco, Q) . Certainly, if we check out all 1' in c whether 
perpendicular or not we would suffer from the heavy computational 
burden. Thus, our main goal is to develop a computational feasible 
algorithm for it . The algorithm is indicated in the appendix and 
our conjectures appear in the appendix also. 
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2.1.2 Family of polynomials with the coefficients of multilinear 
affine function. 

For the family of characteristic polynomials 

P (s, q) = E n i=0 a^qjs 1 V q e Q C R n 

Where a^q) is multilinear affine function for each i, and its zeros 
cr[p(s,Q)] = {z e C : p(z,Q) = 0} 

Our problem is finding necessary and sufficient conditions for 
o tp (s,Q) ] C C*_ , 

That is, p(s, Q) e Hurwitz. 

As we know, the image of Q in the coefficient space, 

a (Q) = {a(q): a(q) = (a c (q), a x (q) , • • •, a n (q)), q e Q} C R n+l 

is not necessarily convex. The stability of a whole family of 
polynomials can not be inferenced from the stability of just the 
set's vertices, edges or boundaries. But consider the box 

A* C R n+1 , 

A* = [A' 0 , A + 0 ] X [A',, A + i] x • • • x [A' n , A + „] , and 

A"; = min qeQ [ a ; (q) ] , A + j = max qcQ [a t (q) ] 

representing polynomials with independently varying coefficients in 
the intervals defined by the minimum and maximum values of the 
multilinear affine coefficient functions, a;(q) , defined on Q. 

For new box A*, we have a family of polynomials 

p ( s , A* ) = { p ( s , A) : p (s, A) =E " i=0 A i s i ,A= (Ao,Aj, • • • ,A„) , ACA*} 

and it becomes interval polynomials. Then we have 

p (s,Q) C p (s, A*) 

While Kharitonov's Theorem guarantees us that the interval 
polynomials in A’ are Hurwitz if and only if Kharitonov's 
polynomials {K n (-), K 12 ( • ) , K 2I ( • ) , K 22 ( • ) } C A* are Hurwitz. 

It is evident that 

p(s,A*) e H p(s,Q) e H 
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Therefore we have the following theorem: 

Theorem 1 . 


p (s , Q) e H if and only if Kharitonov's polynomials 
{*„(•), K 12 (*), K 21 ( • ) , Kjjl-) } C A* are Hurwitz. 
where A* = [A' 0 , A + 0 ] x • • • x [A" n , A + n ] 

and A'j = min q t q a; (g) ; A + j = maXq f g a x (g) . i = 0,1, ...,n. 

2.1.3 Diamond families of polynomials with complex coefficients. 

As we know, in some cases the stability of a whole family of 
polynomials can be inferred from the stability of just the set's 
vertices, edges or boundaries. Here, we consider the case of the 
diamond family of polynomials. We have then a result as following: 

Theorem 1. Let the complex diamond polynomial family be 

W (s) = { p ( s ) : p ( s ) = E n i=0 ( a i + jbi> si 

£"=ot I a r a il + Ibi - £il ) * r} C P" 

Where P n = {p(s) : p(s) = E n i=0 (a^ jb s ) s 1 , a n + jb n 5* 0} 

and a,, £j are the coefficients of the nominal polynomial 

Po ( s ) = E n j=0 < a i + j^s 1 

Then W(s) e H if and only if p,(s), p 2 (s), , p 8 (s) e H 


Where Pi(s) 

= 

Pots) 

+ 

r 

p 2 (s) 

= 

Po(s) 

- 

r 

p 3 (s) 

= 

Pots) 

+ 

jr 

p 4 (s) 

= 

Pots) 

- 

jr 

p 5 (s) 

= 

Pots) 

+ 

rs" 

Pets) 

= 

Pots) 

- 

rn" 

p 7 (s) 

= 

Pots) 

+ 

j rs" 

p 8 (s) 

= 

Pots) 

” 

j rs n 

Theorem 2 . 

Let W(s) = {p(s) 

: p (s) 

= S n i 

=0 

(aj + jbjjs 1 . 

i = 0, 1, • • • ,n 

and E 

Voter 

a i) 

+ (bj-£j) ] 


W' (s) e H if and only if 4 vertices are Hurwitz and 4 line segments 
are Hurwitz. 
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Details please see Appendix. 

2.2 Strictly Positive Real Functions (SPR) for Robust Design. 

We have known that the sets 

SPR{n,(s)} = {d(s) : Refn, ( jco) /d ( jo>) }>0 and 

n;(s) e H for i = 1,2} 

are convex cones, and the convex combination of n,(s) and n 2 (s) 
keeping Hurwitz is only a necessary condition for existence of d(s) 
such that both n,(s)/d(s) and n 2 (s)/d(s) are SPR. Unfortunately, we 
could not know if it is also a sufficient condition. 

For this, we need to describe the set SPR{n(s)} in some way so that 
to prove SPR{nj(s)} n SPR {n 2 (s)} 5 * 0. Even though we have had 5 
definitions for SPR{n(s)}, they would not work for our problem. 
Recently, we find the following new ways to define SPR{n(s)}. 

Definition 1. SPR{n (s) } = {d (s) : (a 0 + joeb 0 , • • • , a n +jocb n ) eH, VoceR} 

where n (s) =E n i=0 a i s i d(s) = E n i=0 b i s i 

Now, we may then answer our question 

SPR{n, (s) }nSPR{n 2 (s) for Xn, (s) + (1-X) n 2 (s) e H. 

Xe [0,1] . 

In fact, consider 

SPR{ Xn, (s) + (l-X)n 2 (s) } = {d(s) : [Xn, (s) + (1-X) n 2 (s) ] +jocd(s) eH VoceR} by 
our definition 1 

Then, we have that Xn, ( s) + (1-X) n 2 (s) +Xjocd (s) + (1-X) joed (s) eH 

Xn, (s) + (1-X) joed (s) + (1-X) n 2 (s) +Xjocd (s) eH 

Xn, (s) + (1-X) joed (s) eH and (1-X) n 2 (s) +Xjocd (s) eH 

Therefore there exists d(s) e SPRfn, (s) }nSPR{n 2 (s) } 

That is, there exists d(s) such that n, (s)/d(s) and n 2 (s)/d(s) are 
SPR. We will develop further results. 
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2 . 3 ACS . 


After reviewing the XTE Subsystem Analysis Results Report, we 
have realized that ACS is a flexible body linear system with time- 
varying, uncertain parameter and with the uncertain disturbance of 
same movements. Therefore, ACS would be in the form 

X (t) = [A (t ) +AA (r (t) ] X (t ) + [B (t) + AB (s (t) ) ]U(t) +D(t) W(t) (1) 

where AA(r(t)) and AB(s(t)) are yielded by the parameters and the 
control U(t) . W(t) is the disturbance of some moments. Thus, our 
framework for analyzing ACS would be the following: 

1. Construct state equations (1) in time domain for ACS. That is, 
we should look for the conditions for A, aA; B, aB; D, W such that 
the control U(t) stabilizes ACS. For the system (1), we usually 
adopt the Quadratic cost functional 

J min = J°°o[Y(t) T y(t) + U ( t ) t RU ( t ) ] dt 

where R is a positive definite 

to find the solution of its corresponding Riccati's equation to 
obtain the robustness conditions for this system. 

2. Determine Stability Margin for ACS. 

That is, to find the best U(t) such that the system ACS 
maintains asymptotic stability in the region of the uncertain 
parameters and the uncertain disturbances. Therefore we need to 
define a radius p(p,u) and find 

Max^u Min p „ {p(p,u)} 

Where u is a control, p is a disturbation 

Thus, H a norm would be considered. 

3. Adjust and revise u(t) according to the reality of ACS. 
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